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Abstract — In this paper we show that the (co)chain complex 
associated with a decomposition of the computational domain, 
commonly called a mesh in computational science and engineer- 
ing, can be represented by a block-bidiagonal matrix that we call 
the Hasse matrix. Moreover, we show that topology-preserving 
mesh refinements, produced by the action of (the simplest) Euler 
operators, can be reduced to multilinear transformations of 
the Hasse matrix representing the complex. Our main result 
is a new representation of the (co)chain complex underlying 
field computations, a representation that provides new insights 
into the transformations induced by local mesh refinements. 
Our approach is based on first principles and is general in 
that it applies to most representational domains that can be 
characterized as cell complexes, without any restrictions on 
their type, dimension, codimension, orientability, manifoldness, 
connectedness. 

Index Terms — Computational geometry, geometric modeling, 
mesh generation, topology, spatial data structures, algorithms, 
sparse matrices, finite element methods. 



I. Introduction 

A. Motivation 

The boundary representation has become the representation 
of choice in many academic and virtually all commercial 
solid modeUng systems. As a consequence, most geometric, 
scientific and engineering applications have to be formulated in 
terms of boundary representations, often leading to nontrivial 
representation conversion problems. Well known examples of 
such problems include Boolean set operations, finite element 
meshing, and subdivision algorithms. 

Formally, all boundary representations are widely recog- 
nized as graph-based data structures [l]-[4] representing one 
of several possible cells complexes [5]-[8]. Space require- 
ments and computational efficiency of such data structures 
have been studied in the literature (see, e.g., [9]). Historically, 
such cell complexes have been restricted to (unions of) two- 
dimensional orientable manifolds, but a number of extensions 
to more general orientable cellular spaces have been proposed 
(see, e.g., [8], [10], [11]). Depending on a particular choice 
of data structures, boundary representations are constructed, 
edited, and updated using a small set of basic operators on 
the graph representation, while preserving and/or updating 
the basic topological invariants of the cell complex. Such 
operators are commonly called Euler operators [3], [10], [12], 

Paper accepted for publication in a Special Issue on Solid & Physical Mod- 
eling of the Transactions of Automation Science and Engineering, tentatively 
scheduled for the July 2009 T-ASE issue. 

A. DiCarlo, F. Milicchio, and A. Paoluzzi are with the University Roma 
Tre (Roma, Italy); V. Shapiro is with the University of Wisconsin-Madison. 



because they enforce the Euler-Poincare formula. All higher- 
level algorithms and applications of boundary representations 
are implemented in terms of such operators. 

This evolutionary development of boundary representations 
also led to several fundamental difficulties: 

• Variety of assumptions about the cell complexes and 
graph representations make standardization difficult. This 
in turns complicates the issues of data exchange and 
transfer, and leads to proliferation of incompatible algo- 
rithms. 

• Boundary representation algorithms are dominated by 
graph searching algorithms (boundary traversals) that 
tend to force serial processing. Nor is it clear how to 
combine such graph representations with multi-resolution 
representations and algorithms. 

« Extending boundary representations to more general cel- 
lular spaces has proved challenging. Despite many pro- 
posals, most commercial systems are still restricted to 
two-dimensional orientable surfaces. 

• Last, but not least, solid modeling has developed into a 
highly specialized discipline that is largely disconnected 
from many standard computational techniques. In par- 
ticular, boundary representations do not appear to be 
directly related to the methods for physical analysis and 
simulation such as finite differences, finite elements, and 
finite volumes. 

In this paper, we claim that all representations of cell com- 
plexes are properly represented by a (co)chain complex [13], 
[14]. It captures all combinatorial relationship of interest in 
solid and physical modeling formally and unambiguously, 
using standard algebraic topological operators of boundary d 
and coboundary S. We show that the (co)chain complex can 
be represented by a sparse block-bidiagonal matrix that we 
call the Hasse matrix. We also show that topology-preserving 
refinements of such cell complexes correspond to simple 
Euler operators and are easily formulated as multi-linear 
transformations of the Hasse matrix. There are at least three 
important consequences of our proposal. First, the proposed 
approach applies to all cell complexes, without restrictions 
on type, dimension, codimension, orientability, manifoldness, 
and so on. Secondly, as we will see in Section Vlll, this 
formalism unifies geometric and physical computations within 
a common formal computational structure. And finally, our for- 
mulation explicitly shows that many geometric computations 
(and computations with boundary representations in particular) 
can be formulated and implemented in terms of standard sparse 
matrix computational techniques, opening a possibility for a 
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wide range of computational breakthroughs and opportunities. 

B. Related Work 

Algebraic-topological properties of boundary representa- 
tions are well understood — see [3], [5], [8], [15] for details. 
Branin [16] and Tonti [17] advocated a unified discrete view 
of all physical theories using concepts from algebraic topology 
and the de Rham cohomology. More recently, this early 
research led to new efforts in developing unified computational 
models and languages for analysis, simulation, and engineer- 
ing design. Notably, Palmer and Shapiro [18] proposed a 
unified computational model of engineering systems that relies 
on concepts from algebraic topology. A number of researchers 
went beyond the use of chains and cochains as general-purpose 
data types, considering that a sound numerical method should 
reflect the algebraic-topological structure of the underlying 
physical theory in a faithful way. Notably, Strang [19] ob- 
served that the FEM encodes a pervasive balance pattern, 
which is at the center of the classification in [17]. Mat- 
tiussi [20] provided interpretations of FEM, FVM, and EDM 
in terms of the topological properties of the corresponding 
field theory. Tonti [21] presented his cell method as a direct 
discrete method, bypassing the underlying continuum model. 
In [22] EDMs that satisfy desired topological properties are 
discussed. In our previous work [23], physical information 
is attached to an adaptive, full-dimensional decomposition of 
the computational domain. Giving pre-eminence to the cells 
of highest dimension allows to generate the geometry and 
to simulate the physics simultaneously. Such a formulation 
removes artificial constraints on the shape of discrete ele- 
ments and unifies commonly unrelated finite methods into 
a single computational framework [24]. Our goal is to graft 
this approach to field modeling onto an already established 
computational framework for geometric modeling with cell 
complexes [25]. This framework has been recently extended to 
provide parallel and progressive generation of very large data- 
sets using streams of continuous approximations of the domain 
with convex cells [26]. The approach also supports progressive 
Boolean operations [27], providing continuous streaming of 
geometrical features and adaptive refinement of their details. 

C. Overview 

This paper is a revised and augmented version of our 
contribution to the 2007 ACM Symposium on Solid and 
Physical Modeling [28], that now explicitly accounts for unit 
(co)chains, cell size and inner product. We also made an 
effort to render the present version more palatable for the 
prospective reader, and as self-contained as possible. All 
technical notions — apart from the most basic ones, listed at the 
beginning of Section II — are explicitly defined and explained 
through examples. 

In Section II we review the basic notions of (co)chain and 
(co)boundary and their representations using matrices and 
the Hasse diagram. We make clear that chains with real 
coefficients serve to render cells measurable; correspondingly, 
cochains represent densities of physical quantities with re- 
spect to cell measures. The duality between boundary and 



coboundary operators and its representation through matrix 
transposition are highlighted. Section III introduces our block- 
matrix representation of a (co)chain complex. In Section IV 
we use algebraic-topological notions to define a minimal set 
of operators as transformations between cell complexes that 
preserve the Euler characteristics. These operators are shown 
to correspond to multilinear transformations of the Hasse 
matrix in Section V. Section VI demonstrates how common 
algorithms for splitting a cell complex may be formulated in 
algebraic-topological terms. Section VII explains how the pro- 
posed representation may unify geometric and physical mod- 
eling in a common computational framework, showing how 
local adjacency information — leading to the discrete Laplace- 
deRham operators — is straightforwardly obtained through sim- 
ple linear algebra, after chains have been identified with (real- 
valued) cochains. The central role played here (and only 
here) by a further independent mathematical structure — the 
inner product between (co)chains — is brought to the fore in 
Sections VII and VIII. Section IX provides some concluding 
remarks. 

II. Background 

We take for granted the elementary notions of simplex, cell, 
orientation, cell complex, face, and refer the reader to [5], [13]. 
(In Eig. I, a2,cT2' and cr| are simplices, is not; the faces 
of ^re al,af,ai, and af.) Apart from 0-cells, we take all 
cells as relatively open. 




Fig. 1. A 2-complex K. Notice that its 2-ceIIs are colierently oriented. 

A. Chains and Cochains 

Let K he a cell complex representing a finite partition of 
a compact set D C E"*, the d-dimensional Euclidean point 
space (or, more generally, D C Ai, a Riemannian manifold 
modeled on E**). By way of illustration, Eig. 1 exhibits a small 
2-complex partitioning a polygonal subset of E^. We call p- 
skeleton Kp C K the subset of oriented p-cells of K, and 
denote with kp the number of p-cells: kp := hence 

#A' = /co + fci + • • • + fed ■ 

Eor the 2-complex K in Eig. 1, e.g., = 6 + 9 + 4. The p- 
skeleton Kp will be ordered by labeling each p-cell ap with 
a positive integer: Kp = {ap , . . . , ap''). In the following, 
the ordinal and/or the dimension of cells will be dropped 
from notation whenever convenient, and the complex K will 
be identified with the tuple of its ordered p-skeletons; K = 
{Kq, . . .,Kd). 
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1) Chains: Let (G, +) be a nontrivial abelian (i.e., com- 
mutative) group. A p-chain of K with coefficients in G is a 
mapping Cp : Kp G such that, for each a G Kp, reversing 
a cell orientation changes the sign of the chain value: 

Cp(-cr) = -Cp(cr) . 

Chain addition is defined by addition of chain values: if 
Cp,Cp are p-chains, then (c^ + Cp){a) ~ Cp{(^) + Cp{a), 
for each a e iiTp. The resulting group is denoted Cp{K; G). 
When clear from the context, the group G is often left implied, 
writing Gp{K). 

Let a be an oriented cell in K and g E G. The elementary 
chain whose value is g on a, —g on —a and on any other 
cell in K is denoted ga. Each chain can then be written in a 
unique way as a (finite) sum of elementary chains: 

Cp = Y.gk<yp. (1) 

k=l 

Chains are often thought of as attaching orientation and 
multiplicity to cells: if coefficients are taken from the smallest 
nontrivial group, i.e., G = {—1,0,1}, then cells can only 
be discarded or selected, possibly inverting their orientation. 
However, we are mainly interested in reaZ-valued chains, and 
strongly rely on the extra structure imparted to the set R of 
real numbers by multiplication. In fact, the ability to add 
(and subtract) chains is not enough: the very idea of mesh 
refinement (and coarsening) requires chains to be also scaled. 
Hence the full linear structure is brought to bear, as argued 
below. 

Real-valued chains attach a signed p-measure to p-cells — 
such as length to 1 -cells, area to 2-cells, volume to 3-cells, 
thus restoring part of the geometrical information left out 
by the purely topological construction of a cell complex. 
Thanks to the multiplicative structure of the real field, all K- 
valued elementary chains on the p-cell (7^ are multiples (or 
submultiples) of the unit chain Up, whose value is 1 on (jp, 
— 1 on — dp and on any other cell in Kp. Therefore, each 
chain Cp can be seen as a linear combination of unit chains, 
and Gp(A';R) as a linear space over M, spanned by the set 
of unit p-chains, which constitutes its standard basis: 

Cp = J2^ku';. (2) 

k=l 

Notice that representation ( 1 ) would formally coincide with 
(2), if each cell (7^ were identified with the corresponding unit 
chain u^. However, this seemingly innocuous move should 
definitely be avoided, being tantamount to assigning one and 
the same size to all cells. On the contrary, we want to tune 
size freely, by identifying each cell with a selected elementary 
chain, i.e., a multiple of the corresponding unit chain: 

fJp = (Up "p , (3) 

the scalar p^ being the size imparted to the cell Cp. The size 
Pp is assumed to be positive (orientation being accounted for 
by Up ) and significantly different from zero:' p^^O . 0-cells 

'in Section II-B. 5 cell sizes will be used as divisors. Cf. also the definidon 
of sgn^ in Section VI-A. 



are systematically given a unit size: /ig = 1 4=> ctq = Uq . 
Our motivation is illustrated in the following straightforward 
example. 

Example 1 (Measuring cell size): Three different partitions 
of the same 1-D interval are arranged in Fig. 2. The upper and 
middle ones are indistinguishable as cell complexes, i.e., on 
mere topological grounds. As chain complexes, however, they 
may differ: their two elementary 1 -chains, if meant to represent 
(cell length) /(interval length), take the values (.5, .5) for the 
upper mesh and (.75, .25) for the middle one. The lower mesh 
is a refinement of the middle one: its four 1 -chains take the 
values (.5, .25, .125, .125). 

ol 

• >m 

1 4 2 
• >• 

al ol ol ol 
• >• 

Fig. 2. Three ID meshes. 

2) Cochains: By definition, the group of p-cochains of K 
with coefficients in G is the group of homomorphisms of 
Gp(/v;G) into G: 

GP{K; G) := Hom{Cp{K; G), G), 

i.e., -fP e CP{K- G) if and only if 7P : Gp{K; G) G and 

k=l k=l 

for all (.gi,...,5fcj eG''". 

If (co)chains are real-valued, then '^^{Xup) = X'^P{up), 
and each cochain 7^ can be seen as a linear combination of 
the unit p-cochains !][,... ,rjj^ , whose value is 1 on a unit 
p-chain and on all others: 

7'' = EAS,r, where ./K) = { J l^^^.;^ (4) 

Therefore, G^(A';]R) can be seen as the dual space of 
Gp{K;M.). It is spanned by the set of unit p-cochains, which 
constitutes the basis dual to the standard basis of Gp{K;M.). 
The evaluation of a real-valued cochain is aptly denoted as a 
duaUty pairing, in order to stress its bilinear property: 

7^(cp) = (7^,cp). 

Real-valued cochains represent densities with respect to the 
measures imparted to cells by real-valued chains, and the 
above duality pairing is a discrete preliminary to integration: 
its value yields the integral of the density over the chain 
Cp, i.e., the content in the support of Cp of the quantity gauged 
by 7P. 

B. Boundary and Coboundary 

1 ) Boundary Operator: The boundary operator 

dp : Gp{K) ^ Gp-iiK) 
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is first defined on simplices (see [13, Ch. 1 §5]). If Op is tiie 
oriented p-simplex spanned by the ordered set of p + 1 points 
(cco, . . . ,Xp), then 

p 

dp Op := ^(-l)''crp_i^ fc , 

fc=0 

where (Tp_i,fe is the fc-th/ace of o-p, i.e., the oriented (p— 1)- 
simplex spanned by (a;o, . . . , a;fe_i, a;fc+i, . . . , Xp). E.g., with 
reference to Fig. 1, we have 82 (t| = — ctj — + af. 

The next step is to extend dp to cells by partitioning them 
into simplices, and exploiting the additivity of the dp operator. 
E.g., the nonsimplicial 2-cell a\ in Fig. 1 may be partitioned 
into (at least) two 2-simplices. Of course, partitioning is not 
unique, but the result does not depend on the partition used, 
since the contributions of all faces internal to the partitioned 
cell cancel out (see [13, Ch. 4 §39]). For the nonsimplicial 
cell in Fig. 1, one gets 82 <j\ = a\+ a'l + (t\+ af. 

The boundary operator is then extended to elementary 
chains, by taking 

dpigcr) ■■=g{dp(j) , 

and finally to all chains by additivity. In conclusion, the 
boundary operator is a homomorphism: 

8p e Hom{Cp{K),Cp.i{K)). 

In particular, it preserves the linear structure real-valued chains 
have been endowed with. In the following (see Section II-B. 5), 
we shall represent 8p by means of the fcp_i x kp matrix [8p\ 
of its components with respect to the standard bases of Cp{K) 
and Cp.iiK). 

2) Coboundary Operator: The coboundary operator 

5P : CP{K) CP+\K) 

acts on p-cochains as the dual of the boundary operator dp+i 
on {p + l)-chains: for all -fP G and Cp+i G Cp+i, 

{SPjP, Cp+i)p+i = (7'', dp+i Cp+i)p . 

For once, we have labelled each duality bracket with the 
dimension of the (co)chains it acts on. Recalling that chain- 
cochain duality means integration, the reader will recognize 
this defining property as the combinatorial archetype of Stokes ' 
theorem. Denoting the dual of an operator by starring its 
symbol, we shall write: 

f)P - 8"^ 
— i^p+i ■ 

By definition, the coboundary operator is a homomorphism 
preserving the linear structure real-valued cochains have been 
endowed with: 

6P e Hom{CP{K),CP+\K)). 

In the following (see Section II-B .5), we shall represent 6^ hy 
means of the fcp+i x kp matrix [S^] of its components with 
respect to the standard bases of Cp{K) and Cp~^^{K), i.e., 
the bases dual to the standard bases of Cp{K) and Cp+i(A'). 
Since we use dual bases, matrices representing dual operators 
are the transpose of each other: for all p = 0, . . . , d— 1 , 



Example 2 (Cochains, Coboundary and Electrostatics): 
Consider once again the 2-complex K in Fig. 1, its support 
[[A'D {i.e., the point set union of all cells in K), and a scalar 
field (f) : |K] M, to be thought of as the electrostatic 
potential. The discrete counterpart of (p is the 0-cochain 0° 
sampling it at the six 0-cells of K. The electric field E is the 
negative gradient of the electrostatic potential: E = —W(f>. 
A measurable quantity is its circulation along any curve 
c : [0, 1] |K], defined as 

C(E,c):= / E(c(s)).c'(s)ds, 
Jo 

where a prime denotes differentiation with respect to the 
curvilinear coordinate. Since E ~ — V(/), it turns out that 

C(E, c) = - V</)|,(,) • c'(s) ds = cf,{c{0)) - 4>{c{l)) . 
Jo 

The analog in ii' of a curve is a 1 -chain, and the discrete 
counterpart of the vector field E is the 1-cochain C ^ sampling 
its circulation along the nine 1 -cells of K, delivered by the 
opposite of the coboundary of </)": 

i.e., is the discrete counterpart of V acting on scalar fields. 
Since E admits a potential, its circulation along the boundary 
of any disk V C |K] is null: 

C(E,aP) ^ I rotE = 0. 
Jv 

Analogously, for any linear combination of the four 2-cells in 
K, i.e., for any 2-chain C2 G K2, we have 

{C\82C2) - {6^C\C2)^0, 

since 5^o S'^ ~ (cf. Section II-C.l). Clearly, is the discrete 
counterpart of the curl operator on vector fields. 

3) Matrix Representation of Chains and Cochains: Since 
p-chains and p-cochains form dual linear spaces (cf. Section II- 
A), they lend themselves to the usual representation of vectors 
and covectors as column and row matrices. The standard 
basis of Cp (cf. Section II- A. 1) and the dual basis of C (cf. 
Section II-A. 2) will be used throughout. 

The components gi, ■ ■ ■ ,gk of a p-chain Cp in (2) may 
be organized into a column matrix Cp ~ [cp]. Analogously, 
the components g^, . . . , g^p of a p-cochain 7*^ in (4) may be 
organized into a row matrix y*' = [7^*]. The duality pairing 
between 7^* and Cp is represented by a matrix product: 

(7ncp) = [7^][cp]. 

4) Incidence Matrices: The intersection between p-cells 
and relatively closed [p + l)-cells may be characterized by 
the p-incidence matrix Bp, whose entries B^^ are defined as 
follows: 

Bp^ =0 if dp n o'p_^,i = % ,a being the closure of a; 
Bp' = ±1 otherwise, with +1 (—1) if the orientation of dp 
is equal (opposite) to that of the corresponding face 

of cr^+i . 

Be it noted that the information stored in Bp is purely 
topological. As is obvious, its transpose Bp describes how 
[p + l)-cells intersect with p-cells. 
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5) Matrix Representation of Boundary and Coboundary 
Operators: We now introduce the notion of measured p- 
incidence, by defining a matrix Mp whose entries M^^ depend 
on both topology (through p-incidence) and measure (through 
cell size): 

It should be no surprise that Mp is exactly the matrix 
representing the boundary operator 9p+i with respect to the 
standard bases of Cp+i and Cp'. 



[dp+i] = Mp . 



(5) 



Consequently (cf. Section II-B. 2), its transpose represents the 
coboundary operator with respect to the dual bases: 



(6) 



of the boundary operators di and 82 and of their duals, the 
coboundary operators = d* and 6^ — 82 is illustrated by 
arrows. 
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->«(4.,0.) 
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Fig. 3. A small 2-complex. Ordered pairs are Cartesian coordinates in E 



If size represents length, area, volume for 1-, 2-, and 3-cells, 
respectively, then measured-incidence matrices have physical 
dimension (length)^^. Hence, [d] and [S] act as first-order 
difference operators. 

Measure may be readily filtered out of Mp, and Bp recov- 
ered, whenever needed: 



p 



sgn(M;^). 



C. Hasse Diagram of a Chain Complex 

Hasse diagrams, named after the German mathematician 
Helmut Hasse (1898-1979), illustrate the cover relation of a 
partial order, and are commonly used for representing partially 
ordered sets — or posets, for short. 

If is a poset, its Hasse diagram is the graph H = {N, E), 
such that for any x,y G N, there exists {x, y) G E \f and only 
if X < y, and there is no z G such that x < z < y. Given 
a d-complex K, let the sets N, E be defined as follows: 

1) N:=KoUKiU---UKd, 

2) E:=EiU---[JEd, with 

3) Ep := {(o-p, (Tp_i) I CTp_i e dap}, 1 <p<d, 

then the graph 'H{K) := (N, E) provides a complete repre- 
sentation of K. Attaching a label from {+1,-1} to the arc 
(x, y) e Ep, denoted sgn(x-, y), suffices to specify the relative 
orientation between the p-cell represented by the node x and 
the (p— l)-cell represented by the node y. A (strictly positive) 
real label is needed to specify the size ratio fip^i{x) / fip{y). 

Given Ti{K) = {N, E), for each x e N define the subset 
of its "children" iV^ C TV: 

1) E'^:={ix,y)\yeN,ix,y)eE}, 

2) N-:={y\yeN,ix,y)eE'^}. 

Then, the boundary of the unit chain u{x) represented by the 
node X may be computed as: 




P = I •crJ •o-f •cr, 

\ .25-1. ^1. / ^1 



p = 



51 



• C7 



Fig. 4. The Hasse diagram of the 2-complex in Fig. 3. Each arc (a;, y) bears 
a non-null real label \{x,y), whose sign specifies the relative orientation 
between x and y, while its modulus |A(3;,y)| specifies the size ratio 



1 ) Chain and Cochain Complexes: An abstract chain com- 
plex C = {Cp, dp) is defined as a sequence of abelian groups 
Cp, paired with homomorphisms dp, p > 1, satisfying the 
relations dp o dp+i = 0, for each p> 1: 



c 



p-1 



Ci ^ Co . 



Analogously, an abstract cochain complex C ~ (C^, S^) is 
a sequence of abeUan groups C, paired with homomorphisms 
5P, p> 1, satisfying the relations o Sp~^ = (p > 1): 



1 ^p^^ 1 



C' ^ C°. 



If a chain complex C consists of linear spaces Cp and linear 
operators dp, then the dual spaces C^ — C* and the dual 
operators 6^ = 9p_^i make up a cochain complex {€*). 

2) (Co-)Chain Maps: Let C = {Cp,dp) andd = {Cp,dp) 
be two chain complexes. A chain map : C ^ C is a p-family 
of homomorphisms 



bp ■ Cp 



Cn 



du{x) 



yeN' 



sgn{x, y)u{y) , 



where u{y) denotes the unit chain represented by the node y. 
The boundary of a general chain follows by linearity. 

Example 3 (Hasse Diagram): Fig. 3 represents a 2-comp- 
lex embedded in E^, comprising two 2-cells, seven 1 -cells 
and six 0-cells. Fig. 4 depicts its Hasse graph. The action 



such that dp o 
commutative: 



c 



_i o dp, i.e., the following diagram is 

<t>p 



dpi 



C 



Cp 
I dp 

Cp-i 



Cochain maps are defined analogously. 
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III. Matrix Representation 

In this section we introduce a block-matrix representation 
of the dual complexes C = iCp,dp) and C* = (0^,6^) 
associated with a decomposition of the computational domain, 
and call it the Hasse matrix. 



A. Block-Matrix Decomposition 

Let K be a d-complex and 'H{K) its Hasse graph. The 
Hasse matrix H(A') is defined by the block structure shown 
in Fig. 5. Its transpose is the Hasse matrix of the dual complex 
K*, whose Hasse graph Ti.{K*) is isomorphic to Ti{K), with 
K* = Kd-p (0 < p < d), where the boundary and cobound- 
ary operators are swapped by duality: 11{K*) = H(A')^. 

Since matrices representing boundary and coboundary oper- 
ators are the transpose of each other, each block in the Hasse 
matrix may be interpreted either way. As shown in Fig. 5, we 
prefer to label blocks with coboundary matrices (and their 
transposes). The reason is that S's acting on cochains are 
the discrete counterpart of the exterior derivative d acting 
on differential forms. For typographic convenience, we shall 
occasionally write dp in lieu of S^. 



ki 
k3 



ko 



Fig. 5. Block structure of the Hasse matrix: the whole scheme holds for d 
odd; for d even, the last block-row should be discai'ded. 

Example 4 (Hasse Matrix of a 2-Complex): Let K be the 
2-complex in Fig. 3, whose Hasse diagram is represented in 



Fig. 4. 


Its Hasse matrix 


is the fci X 
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The left ki x ko — 7 x 6 block equals [S^] 
fcl X /C2 = 7 X 2 block equals [S = [82]. 



, the right 



IV. EuLER Operators 

In solid modeling it is common to refer to Euler operators 
as an independent set of operators that transform a boundary 
representation of a solid into a different one, satisfying the 



Euler-Poincare formula [3], [15]. They may be allowed to 
change the Euler characteristic, whose definition is recalled 
below. 

A. Euler Characteristic 

A well-known invariant associated with a finite d- 
dimensional cell complex K is its Euler characteristic x(if), 
defined as the alternating sum 

X{K) = ko-ki+k2-k3 + --- + (-1)% . 

For polyhedra homeomorphic to the 2-sphere, the Euler char- 
acteristic is ko — ki + k2 ^ V — E + F ^ 2. 

According to the above, the simplest set of independent 
refining (coarsening) operators for a d-space that do not 
change its Euler characteristic has to increase (decrease) by 
one both fcp_i and kp, for 1 < p < d . Therefore, there 
are d elementary refining operators, and the same number of 
elementary coarsening operators. 

1 ) Properties of the Euler Characteristic: The Euler char- 
acteristic is additive under disjoint union: if AI and N are 
disjoint topological spaces {AI O N = (d), then 

X(MU7V) = X(M) + XW. 

More generally, if AI and TV are subspaces of a larger space 
X, then so are their union and intersection, and the Euler 
characteristic obeys the rule: 

X(M UN)^ xW + xiN) - xiM n N). 

Moreover, the Euler characteristic of the product AI x N equals 
the product of the Euler characteristics of M and N: 

xiMxN)^xiM)xiN). 

B. Make and Kill Operators 

The simplest Euler operators that transform a cell complex 
K into another complex K such that xi^) — x{K) add 
(remove) one p-cell and one (p+l)-cell to (from) the complex. 
They will be denoted respectively by (3 and k, from the Greek 
words 'blastos' (construction) and 'klastos' (destruction). 

By definition, the operator (3^ adds a p-cell and a {p + 1)- 
cell to K, thus transforming it into K. The reverse operator 
deletes a p-cell and a (p — l)-cell. In the following, we 
provide detailed information on the construction of refining 
(or make) operators [3^ {p — 0,..., d—1). Coarsening (or kill) 
operators kP (p = 1, . . . ,d), which we will not elaborate on, 
are constructed analogously. 

Let us first discuss how the whole hierarchy of coboundary 
operators transforms under the action of a make operator: 



(g = 0,...,n-l). 



It is easily seen that (3^ affects in a nontrivial way only the 
coboundary operators whose domain and/or codomain change 
under its action, namely: 

1) Sp+i 

2) Sp^i K-> 5p 

3) Spi-^ Sp, 
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Op-1 




as shown by the following commutative diagram: 

(JP+2 _ Qp+2 < (JP+1 i (JP <- 

"^^^ f3P 



Three different computations have to be performed, depend- 
ing on whether the domain only changes (case 1), or the 
codomain only (case 2), or both change (case 3). 

1} Addition of a Column (5jy+i ^ Sp+i): Let the matrix 
[(Sp+i] be mxn; then, the matrix [Sp^i] will be mx (n + 1). 
The column to be added to [Sp+i] pinpoints the {p + 2)-cells 
in K incident on the new cell cTp+i. It is a linear combination 
of the columns of [(5p+i]: 



[6p+i] = [6p+i] 







^ I ri X Ti 






Cn / 



^4 



5i 



"0 



K2=K2 



Fig. 7. First refining step: K = P''\K) ^ {Ko iJ {aD^KiU {af},K2). 



^4 Ji* 



^0^ 



Ki = p\Ki) 




Fig. 8. Second refining step: K S /J^X) = [Kq, KiU{a^} , K2U{al}). 



2) Addition of a Row (Sp-i i~> Sp^i): The row to be added 
to [Sp^i] pinpoints the (p— 1) -cells in K incident on the new 
cell (7p . It is a linear combination of the rows of 

/ I m X m 



[Sp^i] = -R[6p_,] 



3) Addition of a Column and a Row (5p t-^ Sp): One row of 
[Sp] (one chain in Cp) is substituted by two rows (two chains 
in Cp), whose components on the new cell ap sum up to zero. 
The matrix [Sp] is obtained as the sum 

^ 3 
1=1 

where the first term (i — 1) provides the contribution of the 
split cell (Jp+i, the second (i = 2) the contribution of the new 
cell CTp+i, and the third (i — 3) the contribution of all other 
{p + 1) -cells in K. 

C. Examples 

A minimal 2-complex K is represented in Fig. 6. Cell sizes 
are assigned as follows: /i} = ^J2,^ = /if = 1- ,/^2 = -S- 
Fig. 8 shows its refinement K, obtained by applying first the 
operator ff^ to halve the 1-cell a\, then the operator to 
halve the 2-cell a\. Fig. 7 shows the intermediate stage {K). 



On 



Fig. 6. Coarse complex K = {Kq, Ki, K2)- 




Example 5 ( Coboundary Sq : (A') (K) ): Domain 

and codomain have both dimension 3. From Fig. 6 it is seen 
that the matrix representation of Sq is 

-1/V2 1/^2' 

[So] = \ -1. 1. 
0-1. 1. 



C^{K)): In this 



Example 6 (Coboundary Si : C^{K) 
case, ki = 3 and k2 ~ 1, so that 

[<5i] = (-2V2 2. 2.). 

Example 7 (Coboundary Sq : C"(/v) C^{K)): We have 
now fco = fci = 3+1. In Figs. 7,8 the new cells are highlighted 
in boldface. Since both domain and codomain of Sq are larger 
than those of Sq, the new coboundary operator is given by the 
sum of three contributions (see Section IV-B. 3): 



[So] = ( Si S2 S3 




where 



Si 



/ 


1 









/ 













































, S2 = 













V 








) 




V 


1 





/ 













1 











1 


V 





/ 



Matrices Si, S2 extract the row of [Sq] that corresponds to the 
1-cell al to be split (recall that a row of [So] equals a column 
of [di]): Si associates that row to al, while S2 associates it 
to the added cell af. Matrix S3 selects the remaining rows 
of [So] and keeps them unchanged. Matrices Si act on [So] as 
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follows: 



same operations as in Example 7, we get (see Fig. 8): 



Si [So] 





V / 



S3 [So] 



, S2[So] 










Each column of a matrix corresponds to a 1-cell in K. 
Post-multiplication yields: 



Si [<5o] Ti = 



/ 


-V2 








V2 \ 


































V 














) 


/ 














































V 








72 


-72 J 


/ 











' ) 






-1. 


1. 












- 


-1. 


1. 







V 











) 





S2 [^o] T2 



S3 [<5o]T3 = 



Summing up, we obtain (see Fig. 7): 



[So] = 



Example 8 (Coboundary Si : C^{K)^^ C'^K)) 
case, fci = 3 + 1 and ^2 = 1 . Hence, [Si] = [Si]C 



f-V2 








72 \ 


-1. 


1. 











- 1. 


1. 





V 





V2 


-72/ 



[Si 



.5 .5 
10 
10 



-V2 2. 2. -V2). 



^ Example 9 (Coboundary So ■ C"{K) C^{K)): Here 
fco = fco = 4, fci = fci + 1 = 5, and we have: [So] = R [So] 



( 



\ 





1 


72 



/-72 
-1. 




V 





1. 

-1. 



V2 






1. 

72 





72 




-72 
72 



-72 

-1. 









1. 

72 



72 




-72 y 



Example 10 ( Coboundary Si : 

fci = fci + 1 = 5 and k2 = ^2 



^- 1 = 2. 



-> C^iK)): Here 
By performing the 



Ti 

[Si] = ( Si S2 S3 ) [,5i] ( T2 

T3 



-2^2 4. 2^2 \ 
4. -2^2 -2^2 ]■ 



V. Hasse Transformations 

Let K be a d-complex and J1{K) its nxm Hasse matrix. 
Notice that x(if) = m — n so that, while both m and n in- 
crease under topology-preserving refinements, their difference 
does not. We introduce the Hasse transformation 



H(A') ^ H(/v) 

induced by the make operator /3 that splits the (p 
(Tp_i_j into two cells, namely 



l)-cell 



o-p+i and 
and adds a new p-cell, namely 

~fcp+i 



>+l ' 



Let us assume — to be specific, but without loss of 
generality — that d — Z. In this case, H(X) is comprised of 
two diagonal blocks, [So\ and [^2], and one upper-diagonal 
block, (see Example 4 in Section III). 

Example 11 {Make Operators in 3D): Three different al- 
gorithmic patterns arise, depending on the order of /3^: 

\ 











R[<^i]* 






(Si S2 S3)[5o] 


(1) 




/3°(HW) = 






In this 








[^2] 




I 






r(H(A')) = 



VI. Hyperplane Splitting 

In this section we discuss the SPLIT subdivision algorithm 
introduced in [29], rephrasing it in terms of the algebraic 
machinery developed in the previous sections. Subdivision 
algorithms allow to represent a smooth surface via the gener- 
ation of coarser piecewise-linear meshes. They are the result 
of an iterative process of subdividing each polygonal face into 
smaller faces that better approximate the smooth surface. Sub- 
division is also used in combination with multigrid algorithms. 



9 



for generation of finer 3D domain meshes in computational 
science. 

The algorithm [29] works efficiently on a single d-cell of a 
d-complex, by locally splitting it with a hyperplane. The SPLIT 
algorithm performs first a numerical computation classifying 
all the vertices with respect to a given hyperplane, then carries 
out symbolic manipulations on higher dimensional cells, in 
order to obtain the resulting topology. Our algebraic formula- 
tion is general and easy to implement using standard packages 
for sparse-matrix computation [30]. The SPLIT algorithm is a 
useful tool for refining cell complexes, providing the ability to 
compute Boolean operations when combined with BSP trees 
in a progressive way [27]. The SPLIT algorithm is also useful 
to approximate continuous maps between cell complexes. 

Subdivisions of cell complexes are formally defined as 
follows [13]: K is a subdivision of K if 

1) for each a e K, there exists a G K such that cf C cr; 

2) for each a E K, there exists a finite subset {ai} C K, 
such that (T = Ui CTi. 

The SPLIT algorithm — as detailed in the following Sec- 
tion VI- A — satisfies Property 1 by construction. Property 2 is 
also satisfied, since every cell in K is mapped into the union 
of at most two cells (5^ and 5+) in K. 

A. The SPLIT Algorithm 

Let us first introduce two auxiliary operators, useful for the 
matrix formulation of the SPLIT algorithm. 

Definition 1 (Sign): The function sgn^ : {-1, 0, 1}'* 

returns the matrix listing the signs of the elements Vi of a 
d-tuple v= (vi), accounting for a numerical tolerance £ > 0: 

r -1, V, < -e 
(sgn^v) . = I , -e < Vi < £ 
[ 1, V, > £ 

Definition 2 (Absolute value): The function abs acts on a 
real- valued matrix M = [mij], returning the matrix of the 
absolute values of its elements: 

abs (my ) = (|my|) . 
Let the splitting hyperplane 

t) - {(xi, ... eE-^l ^/i^x, = 5} 

i 

be identified with the linear {i.e., affine homogeneous) form 
_^ characterized by the row matrix 

h.= {h^ ... h'^ -b). 

Let V be the column matrix listing the homogeneous coordi- 
nates of the 0-cell ctq : 

Clearly, ao belongs to the above (below) subspace ()+ ( f)^) if 
and only if hx < : computing the sign of the scalar product 
hx solves the point location problem. After introducing the 
matrix 

X := ( xi X2 • • • Xfeo ) , 



collecting the homogeneous coordinates of all the 0-cells in 
A'o, these are classified with respect to the f) splitting hyper- 
plane by the classification function : A'o — > { — 1,0,1} 
represented by the row matrix 

cO :=sgnJhX). 

The SPLIT algorithm proceeds hierarchically from 0-cells 
up to d-cells by (a) classifying the cells with respect to 
the splitting hyperplane, and (b) updating the cell complex 
accordingly, including the new elements in the skeletons of 
all orders. The algorithm is sketched in Fig. 9. The inequality 
in step 5 characterizes the p-cells that intersect the splitting 
hyperplane t). 



algorithm split (input: K, X, h; output: K, X); 

1) p:=0 

2) Classify 0-cells: cO:=sgn^(hX) 

3) p:=p+l 

4) Classify p-cells: 
cP:= (absB^_i)cP-i 

aP:= (absB^_i)abscP-i 

5) foreach |cf | 7^ do: Update the cell complex: 

split the i-th p-cell: K := pP-^{K); 
set the new element value: c?."^ :=0 

hp— 1 

6) Re-classify the p-cells of the updated cell complex: 
cP:=sgn,((absB^„i)cP-i) 

7) if p < d then goto step 3, else stop . 

Fig. 9. The SPLIT algorithm, implemented by using classification functions 
and the p-incidence matrices introduced in Section II-B.4. 



B. Splitting Example 

Let us go back to the 2-complex considered in Section IV- 

C, and refine it with the splitting hyperplane specified in 
Fig. 10(a). The reader should recall Figs. 6-8 and refer to them 
to locate cells by name. 




Fig. 10. (a) The splitting hyperplane 1) ; (b) the classification of 0-cells. 

The SPLIT algorithm is initialized by setting p = and 
classifying 0-cells: 

c" =sgn,(h(xi X2 X3)) = (-1 1), 

as shown in Fig. 10(b). Then, p is increased to 1 and 1-cells 
are classified by computing: 

= (absB^)c" = (0 -1 1), 
= (abs B^) abs c" = (2 1 l). 
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The result is illustrated in Fig. 11 : we see that al should be 
split, since |cj| 7^ aj. The operator adds a new 0-cell 
(classified to 0) and a new 1-cell (see Fig. 12). The two 1-cells 
resulting from the split are reclassified, as shown in Fig. 12. 
Then, p is increased to 2 and 2-cells are classified: 

t 



= (absB^)ci 



1111-1-1 1 1 



0, 



= (abs ) abs = ( 1 1 1 1 ) ( 1 1 1 1 ) = 4 , 

(see Fig. 13a). Hence, ctj gets split, the splitting being exe- 
cuted by the (3^ operator that creates one 1-cell and one 2-cell 
as shown in Fig. 13. Finally, 2-cells are reclassified and the 
algorithm terminated, since p ^ d . The result is illustrated in 
Fig. 13(c). 



+1 



Fig. 11. The classification functions and used to detect 1-cells 
intersecting the splitting hyperplane. 





+1 



+1 



Fig. 12. The updated cell complex, with 1-cells reclassified. 




-1, 



+ 1 




-1 



Fig. 13. (a) Classification of 2-cells; (b) the classification function on the 
refined 1-skeleton; and (c) the refined 2-skeleton with the classification of 
2-cells. 



C. Subdivision of a Complex 

Since SPLIT is a subdivision generator, the process can be 
iterated any number of times: for any n G N, (SPLIT„ o • • • o 
SPLlTi)i^ is a subdivision of the complex K, which we shall 
call an n-iterated subdivision and denote briefly by SPLIT"X. 

From the ^«ife approximation theorem [13] the following 
result follows. Let ip : \K\ |L] be any continuous map 
from the support of a finite complex K to the support of a 
complex L (recall from Example 2 that [[K]] is the point set 
union of all cells in K). Then, for an n-iterated subdivision of 
K with n big enough, there exists a map ip : SPLIT" A' L 
approximating if. 

Another important property of the SPLIT subdivision is 
guaranteed by the algebraic subdivision theorem [13]. The 



splitting induces a unique chain map (as defined in Section II- 
C.2) C : K ^ K := SPLIT K. Consequently, boundaries in 
the refined cell complex K may be evaluated by applying the 
chain map ( to boundaries in the coarse cell complex K. 

VII. Adjacency and Laplace-deRham Operators 

In graph theory, the adjacency matrix of vertices is one of 
the many representations of a graph, i.e., a 1 -complex K ^ 
{Kq,Ki). The well-known relation between incidence and 
adjacency matrices of a graph can be extended to incidence 
and adjacency matrices of all orders p < d of a general d- 
complex, for any d S N . 

Definition 3: Let Mp be the measured p-incidence matrix 
introduced in Section II-B. 5. Post- or pre-multiply Mp by its 
transpose, obtaining respectively: 



A+:=MpM^ 



Mi M„ 



(7) 



The (symmetric) matrix A+ is, by definition, the adjacency 
matrix between p-chains through {p + l)-chains; analogously, 
^p+i is the adjacency matrix between {p+ l)-chains through 
p-chains. 

As established by (5) and (6), Mp represents the boundary 
operator dp+i with respect to the standard bases of Cp+i 
and Cp, and its transpose Mp the coboundary operator 
with respect to the dual bases of C*^ and C^^^. Therefore, 
their products in (7), while legitimate as matrix operations, 
cannot possibly represent products of boundary and cobound- 
ary operators, unless chains are identified with cochains. This 
identification is performed by introducing a sequence of linear 
isomorphisms Gp (0 < p < d) between chain spaces and their 
dual cochain spaces: 



fjp+1 



CP ^ CP- 



— 1 



6° 



%+2 



Gp 



Gp—i 
9. 



p+l > Op > Cp-i 



Cp-1 



9i 



The (nondegenerate) bilinear form induced by Gp on 
Cp will be denoted by the same symbol: for all Cp,c!^ ^ Cp , 



Gp{cp, Cp) .— {GpCp, Cp) 



(8) 



The inner product (8) is assumed to be symmetric: for all 
Cp, Cp G Cp , Gp{c'p, Cp) = Gp{cp, Cp) G* = Gp . 

The adjoint (or transpose) of the boundary operator 
dp+i is characterized by the property: 

Gp(cp, 9p+iCp+i) ~ Gp+i(^dJ^ 1 Cp , Cp^ 1 j , 

to be satisfied for all Cp G Cp and Cp+i G Cp+i. It is easily 
checked that 



p+i 



(Gp+i) ^ dp_^iGp ~ (Cp+i) ^S^Gp. 



The isomorphism Gp is represented by the (symmetric) 
Gram matrix Gp, whose entries C^^ are the components of 
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the Gp-images of the elements of the standard basis of Cp in 
the dual basis of C^: 

GpUp ~ ^ Gp^rjJ . 

Consequently (recall (6)), 

J = G-l,mlGp . 

Therefore, the adjacency matrices introduced in (7) represent 
the (symmetric) adjacency operators 

only if the inner product hierarchy is trivial: for all < p < d , 



(9) 



The discrete Laplace-deRham operators — a fundamental 
and ubiquitous ingredient of physical modeling — are defined 
as sums of duals of adjacency operators: for all < p < d, 



Ap:= {dp+idj 



P+i 



dJdpY = sJSp 



5p-i 



it being intended that do,dd (and hence ,dj) are null. Of 
course, the straightforward representation 



- A: 



is only valid if the inner product hierarchy is trivial. 

Example 12 (Adjacency Matrices in 3D): Let us consider 
the chain and cochain complexes associated with the 3- 
complex K whose topology is given in Fig. 14. As far as 
measure and inner product structures are concerned, we make 
the trivial choices: fi^ = 1., GpU^ = r]f for all < p < d and 
1 < i < fcp . As the reader will notice, these choices — though 
legitimate — do not mimic, not even roughly, the Euclidean 
properties of [[X]] suggested by Fig. 14. 

Boundary and coboundary operators are represented by the 
matrices: 



/-I 

-1 



[^o] = [di 



1 \ 
10 



1 



0-1 
-1 1 
0-1 






1 



1 

-1 1 



0-101 

y -1 1 J 



/-I 1 1 \ 



-1 1 



1 







1 





1 



-1 



1-110 
10 0-11 








-1 -1 
y 1 1 





1 

-1 / 



[^2] = [^3]* = 



11110 
-1 1 1 1 



A straightforward computation 
Laplace-deRham matrices: 

/ 



yields the adjacency and 



[Ao] - [So] 



[So] [di 



1 1 



1 -1 
-1 




[^2] [Si 








[^3] [-52 



[Si] [92 





V 



3 
1 
1 
1 

V 
1 1 



\ 



1 
1 






-1 1 



1 





0-1 

1 

V 0-1 1 



-1 -1 
-1 1 



1 1 

1 -1 



/ 1 1 1 1 
1111 
1 1 1 

1 1 2 

-1 
-1 

V -1 



-1 -1 



3 

1 3 -1 
1 -1 3 
1 -1 -1 
1 
1 
-1 



-1 
-1 

4 
-1 
-1 

1 

-1 



-1 
2 

-1 
-1 1 

1 



-1 
-1 
-1 
4 
-1 

-1 






1 -1 





-1 
1 
1 

3 
-1 


1 -1 

-1 -1 -1 





-1 



1 

-1 

1 -1 3 
1 -1 -1 



[A3] = [<52] [83] 



4 -1 
-1 4 



\ 

-1 
-1 
-1 
3 



\ 


-1 
1 
1 


1 
1 

2 J 

\ 





-1 



-1 

-1 -1 

2 -1 

-1 2 ; 



\ 




-1 
1 
1 

1 / 

\ 



-1 
1 

-1 
-1 

3 



The computation of [Ai] and [A2] is left to the reader, who 
is invited to pay attention to the subtle cancellation of terms. 

VIII. Geometry & Physics Modeling 

The real-valued (co)chain-complex formalism and the as- 
sociated Hasse-matrix representation fit physical modeling in 
a natural and straightforward way. Chains may be tuned to 
endow cells with the physical measures of interest: length, 
area, volume, but also mass, charge, and so on. Cochains, 
on the other side, represent densities of all physical quan- 
tities associated with cells through integration with respect 
to those measures. The coboundary operator is behind the 
basic structural laws (balance and compatibility) involving 
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Fig. 14. A sample 3-complex K ^ {Ko,Ki, 7^2, ^3) ■ 



K2 



physically meaningful cochains [17], [18], [31]. It is also well 
known that p-cochains are the discrete analogue of differential 
p-forms [32], [33]. Correspondingly, the cochain complex 
introduced in Section II-C.l is the discrete precursor of the de 
Rham complex [34]-[36], naturally represented by the Hasse 
matrix (or its transpose). This view on physical modeling has 
been progressively advocated [22], [36], [37] as a way to 
increase numerical stability and accuracy of various numerical 
methods. 

The Hasse matrix H(i^), introduced in Section III- A, pro- 
vides a compact representation of boundary d and coboundary 
6 operators, acting respectively on chains or cochains defined 
on K. This representation depends only on: (i) the topology 
of the underlying cell complex, and (ii) the measure imparted 
to cells by chains. 

A supplementary inner-product structure establishes a 
bridge between the chain and cochain complexes, as needed 
for introducing the adjacency and Laplace-deRham operators 
(cf. Section VII). The underlying topology stays untouched. 
Cell measures may — but need not — ^be related with the inner- 
product structure. The trivial choice (9) we made in Exam- 
ple 12 for the sake of argument, while expedient to use on any 
given K, is totally unrelated — in general — to the geometric 
properties of [if ]] relevant to the physics under consideration. 
Dealing properly with the identification between chains and 
cochains is essential for importing into the model the relevant, 
physics-based inner-product structure. This issue is also basic 
to gain the possibility of a meaningful information transfer 
from a cell complex K to any of its refinements K (and vice 
versa), and to establish a decent notion of convergence for 
refinement sequences. 

A deeper discussion of this point — i.e., how to construct 
physically meaningful inner products between chains — is be- 
yond the scope of the present paper. Here, we just stress that 
the very same data structures and algorithms can be used con- 
currently for solid modeling and physics-based simulations. 

From our vantage point, boundary representations and finite 
element meshes appear as two different aspects of the same 
Hasse representation. Furthermore, there is no fundamental 
distinction between different types of approximation methods. 



as we showed in [23], [24] for linear problems. Within our 
framework, the SPLIT algorithm described in Section VI-A 
becomes a powerful tool for progressive refinement not only 
of shapes, but also of the representation of fields living on 
those shapes. 

IX. Conclusion 

Historically, the development of boundary representation 
schemes in solid modeling was driven by limited compu- 
tational resources, and the usual space-time trade-offs [9]. 
Boundary representation was the solution of choice in order 
to (a) save memory, when RAM was small and expensive, and 
(b) spare disk access times, by answering topological queries 
efficiently. 

Contrary to what might appear at first sight, the present 
approach does not imply higher theoretical complexity, since 
the number of non-zero elements in the Hasse matrix H(A') 
is essentially of the same order as the number of adja- 
cency pointers in a typical graph-based representation of the 
cell complex K. Furthermore, the Hasse matrix serves as a 
unifying standard for all boundary representations, different 
graph structures corresponding to different methods [30] for 
encoding the nontrivial information contained in the sparse 
matrix H(/v). 

It should also be noted that chain complexes are a standard 
tool for representing and analyzing topological properties of 
arbitrary cellular spaces. It follows that the proposed Hasse 
matrix and Hasse transformations may codify much more gen- 
eral models, without restrictions on orientability, (co)dimen- 
sion, manifoldness, connectivity, homology, and so on. The 
resulting framework, centered on a matrix representation of 
the domain of interest, unifies several geometric and physical 
finite formulations, and supports local progressive refinement 
and coarsening. 

This approach is strongly motivated by the applications 
to be developed within the next generation of computational 
sciences, which demand large-scale simulation models of field 
problems where geometric and physical properties have to be 
generated, detailed, and refined simultaneously and progres- 
sively. 
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